Fluorescence antibunching microscopy 
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Breaking the diffraction limit in microscopy by utilizing quantum properties of light has been 
the goal of intense research in the recent years. We propose a quantum superresolution technique 
based on non-classical emission statistics of fluorescent markers, routinely used as contrast labels 
for bio-imaging. The technique can be readily implemented using standard fluorescence microscopy 
equipment. 



o : 

<N ■ 
Oh 



CZ3 . 
O . 

• i-H , 

q : 
o ■ 

C/3 ■ 
>Y 

^ : 



(N 

> 

o 
in 



X 



Increasing the imaging resolution in optical microscopy 
can potentially benefit many fields of research, including 
life sciences. In classical linear optics, diffraction im- 
poses a limit on the resolution of far-field microscopy. In 
the last two decades, a number of techniques have been 
developed that break this limit by making use of non- 
linear optical processes [l~3j. The diffraction limit can 
also be effectively overcome by utilizing contrast labels 
exhibiting strong variations in brightness, either induced 
by photoactivation0-ll] or intrinsic @jlH> in effect allow- 
ing for localization of individual fluorophores while the 
rest are dark. 

Quantum optics offers another promising pathway to 
superresolution imaging. Quantum optical methods have 
been shown to dramatically increase the resolution in 
interferometric measurements and allowed for imag- 
ing sensitivity enhancement beyond shot noise [Ioj|. At 
the same time, superresolution via quantum imaging has 
not yet been demonstrated experimentally. The theo- 
retical research aiming at achieving sub-diffraction lim- 
ited quantum imaging has mainly focused on the scenario 
wherein an absorptive object is illuminated with a non- 
classical light beam. Superresolution is then attained 
via two main routes. One requires an object consist- 
ing of or stained with a multi-photon absorbing mate- 
rial. In this scheme, the diffraction limit is overcome 
by utilizing the high spatial frequency quantum inter- 
ference patterns, similarly to quantum superresolution 
lithography [11]. Although multi-photon interference pat- 
terns of high order have been observed experimentally 
using coincidence detection [l^, 13 1, the lack of a low 



light level multi-photon absorber makes this approach 
currently not feasible. The other pathway to achiev- 
ing sub-diffraction resolution is probing an object with 
a beam of light exhibiting position (or momentum) en- 
tanglement. The high resolution images can then be ob- 
tained simply by coincidence detection [l4|, EH . Although 
momentum entangled light can be produced by sponta- 
neous parametric down conversion 16, 17], the resolution 
in this case is limited by the diffraction limit at the pump- 
ing wavelength. Increasing the resolution further requires 
a bright highly entangled light source, which is yet to be 
developed. 

The above superresolution schemes exploit quantum 
properties of the illuminating light and thus require non- 



classical light sources. The alternative is, to image an 
object naturally emitting non-classical light, such as cas- 
caded emission from a three-level system [18j or resonant 
fluorescence (l9j. The multi-photon interference patterns 
observable in the far field can be used for superresolved 
imaging of the emitters. While this approach can be 
feasible for imaging trapped ions, its reliance on frag- 
ile quantum effects makes it impractical for bio-imaging 
conditions. 

In this paper, we consider a different property of fluo- 
rescence emitters: photon antibunching[20|, arising from 
the tendency of fluorophores to emit photons one by one 
rather than in bursts. Antibunching is a distinctively 
quantum phenomenon, implying reduced quantum fluc- 
tuations (squeezing) of light [2l| and sub-Poissonian pho- 
On the other hand, it is a very ro- 
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bust effect, exhibited by various fluorophores at room 
temperature|23-25[. We study the non-classical photon 
statistics of fluorescence in connection to fluorescence mi- 
croscopy and show that it can be used for superresolved 
microscopic imaging under realistic conditions. 

For simplicity, we focus on the case of pulsed excita- 
tion, with the pulse duration much shorter, and the in- 
terval between pulses much longer than the fluorescence 
lifetime. Upon excitation, a single fluorophore emits at 
most one photon, with a probability p. Such behavior 
is profoundly different from the Poissonian statistics of 
classical light. In particular, the variance of the number 
of fluorescent photons emitted in a series of M excitation 
cycles is V — Mp(l —p), while the mean photon number 
is (N) = Alp. A classical light source with the same pho- 
ton flux would yield a variance Vp = (N). The variance 
of the fluorescent photon number is thus reduced by a 
factor of (1 — p) with respect to the classical shot noise. 

The nonclassical statistics of fluorescent light can be 
used to produce superresolved images of fluorophores. 
Consider a fluorescent emitter, imaged by a microscopic 
imaging system onto a pixelated (or scanning) detector 
with a photon number resolving capability. At every de- 
tector position x, the probability P to detect a photon 
emitted by the fluorophore is given by 



P{x) = pQSh(x — xq), 



(1) 



where Q is the quantum efficiency of the detector, S de- 
notes optical collection efficiency, h is the point spread 
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function of the imaging system and xq describes the co- 
ordinates of the fluorophore image. The photon num- 
ber variance becomes a function of detector position: 
V(x) = MP(x)(l-P(x)). For a set of several flu- 
orophores, since the emission events in different fluo- 
rophores are uncorrelated, the variance of the total pho- 
ton number is given by the sum of variances (??) for 
every emitter: 

V(x)=M'£P a (x)(l-P a (x))= (2) 

a 

where p a and Xq are the emission probability and the 
image position of the fluorophore a. 

The difference between eq.© and the classical shot 
noise variance at the same mean flux quantifies the degree 
of antibunching of fluorescent light and can therefore be 
called the antibunching signal: 

A(x) ee V(x) - (N(x)) = ~MQ 2 S 2 Y,Plh 2 {x - x a ). 

a 

(3) 

This signal corresponds to an effective point spread func- 
tion h,A{x) = h 2 (x). In Fourier domain, Ka spans the 
interval of frequencies twice as large as that of h. The 
antibunching microscopy thus enables imaging with up 
to double resolution, similarly to the resolution improve- 
ment attainable with 2-photon microscopy. 

The mechanism of the antibunching microscopy is il- 
lustrated in Fig. [TJ for the case of two identical emitters. 
The emitters are not resolved in the fluorescent signal, 
while two separate peaks are visible in the antibunching 
signal. 

It is instructive to consider photon statistics in the 
limit of small photon flux. Let Fi(x) and F 2 (x) be 
the probabilities of detecting exactly one and two pho- 
tons, respectively, at a given detector position. In terms 
of these probabilities, the average photon number and 
the variance become (N) — M (Fi(x) + 2F 2 (x)) and 
V(x) = M (Fx + AF 2 - Fl) . The antibunching signal 
© then takes the form of A{x) = M (2F 2 (x) - (x)). 
Since this expression vanishes for Poissonian statistics, 
the antibunching signal can be regarded as a measure of 
the lack of two photon coincidence events with respect to 
classical light. 

This observation elucidates the mechanism of the reso- 
lution increase shown in Fig.[TJ in this example, a coinci- 
dence event involves detection of one photon from each of 
the fluorophores. The probability of a coincidence event 
therefore has a sharp maximum positioned between the 
two emitters. This is in contrast to the case of two clas- 
sical emitters, for which a pair of photons could as well 
originate from a single emitter, making the maximum less 
sharp. 

The antibunching signal ([3]) is determined by the op- 
tical signal autocorrelation at a given detector position. 
Fluorescence antibunching is also manifest in the cross 




FIG. 1: Antibunching imaging of two fluorescent emitters. 
The vertical lines denote the locations of two fluorescent emit- 
ters. The blue line is the regular fluorescent signal, the cyan 
line represents the variance given by eq.([2|, the green line 
is the antibunching signal The solid red line shows the 
probability F2 of a two-photon coincidence event with a sharp 
peak in the center. The dashed red line shows the same prob- 
ability for two classical emitters, featuring a wider peak. The 
plots are calculated for Pi = P2 ~ 0.4, with a gaussian point 
spread function with an RMS width of a = 1.06. 

correlation between the photon numbers N\, N 2 de- 
tected in a pair of proximate detectors. Similarly to 
variance, the covariance of the two signals V x (xi,x 2 ) — 
(N±N 2 ) — (Ni) (N 2 ) is a sum of individual fluorophore 
contributions. For classical light, the signals observed in 
the detectors are uncorrelated and therefore covariance 
vanishes. In contrast, an individual fluorophore produces 
only one photon at a time, which can be detected in only 
one of the detectors, leading to NiN 2 = 0. The covari- 
ance antibunching signal can thus be defined as: 

A*(x 1)Xi ) = "£v*(x u x 2 ) = -M^P«P( 2 ) (4) 

where x\,x 2 are the detector positions, P& and P^ 2 ' are 
the detection probabilities for a given fluorophore at the 
two detectors. The antibunching signal (j4]) contains a 
product of two optical point spread functions. As we 
show below, it can be used to produce higher resolution 
images along with the autocorrelation antibunching sig- 
nal ©. 

We numerically tested the resolution improvement in 
antibunching microscopy by performing a Monte Carlo 
simulation of the antibunching imaging process using a 
pixelated detector array. For efficient signal utilization, 
we used both the autocorrelation antibunching signal of 
Eq. ([3]) and the cross correlation data (J4j) to form the im- 
age. The autocorrelation antibunching signal at a given 
pixel, numbered by a two-dimensional index j, is given 
by 

A j =J2(ti) 2 = V j -(N j ), (5) 

a 



3 



where (Nj) is the mean number of photons, and Vj is the 
variance of the photon number. The cross correlation 
contribution to the second order antibunching signal was 
calculated as a weighted sum of the cross correlations of 
the pairs of pixels jiS, centered at the pixel j. The sim- 
ulation was carried out as follows: a stack of frames was 
generated; in each frame, for every fluorophore it was de- 
cided at random whether it emits a photon (with a proba- 
bility p = 0.5). If a photon was emitted, it was randomly 
positioned with a probability density corresponding to a 
Gaussian point spread function. The resulting stack of 
frames was used to compute the second order antibunch- 
ing signal according to the following formula: 

A 3 = Y,W(6) ({Nj+sNj-s) - (N j+S ) (Nj-s)) - (Nj), 
s 

(6) 

where Nj is the number of photons detected in the pixel 
j, 5 is a summation index labeling the pixel pairs, W{5) 
is the weight assigned to the pixel pairs (VF(0) = 1), 
and the angular brackets denote averaging over the set 
of frames. The above analysis does not fully utilize the 
cross-correlation information: indeed, only a half of all 
pixel pairs are centered in a certain pixel. A pair of, 
for example, two adjacent pixels has its center between 
the two pixels. It is therefore possible to compute the 
antibunching signal in 'virtual pixels' in between adjacent 
pixels, i.e. with at least one of the two components of 
j being half-integer 26]. The effective number of pixels 
in each direction is thus doubled, increasing the total 
amount of pixels by a factor of four. The signal in the 
'virtual' pixels was calculated using Eq. (J6]), with the 
last term omitted and with one or both components of 
the summation index 5 assuming half-integer values (so 
that i ± 8 are integer). The results of the simulation 
shown in Fig. [2] demonstrate a significant improvement 
of resolution: the individual emitters, which cannot be 
discerned by regular imaging, are clearly resolved in the 
antibunching image. 

The antibunching signal, defined above in terms of sec- 
ond order momenta, serves as a measure of the lack of 
two-photon coincedence events, hence it can be called 
the second order antibunching signal. The n-th order 
antibunching signal A n (x), quantifying the lack of n- 
photon events, can be defined in terms of the irreducible 
parts of the n-th order momenta known as cumulants |27j]. 
The defining property of cumulants is that they are ad- 
ditive for independent random variables, which allows 
one to express the cumulants of the observed signal as 
a sum of the individual fluorophore contributions. For 
a single fluorophore, the n-th order cumulant is given 
by C n = MP(1 - P) ... (1 - (n - 1)P). The cumulants 
C n (x) of the total signal are therefore given by 

C n {x) = Mj2 p c*{x) (1 - P*(x)) - (1 - (« - l)Pa{x)) . 

a 

(7) 
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FIG. 2: A simulation of the second order antibunching imag- 
ing: resolving individual emitters, (a): seven emitters are ar- 
ranged in a ring. The scalebar shows full width at half max- 
imum (FWHM) of the optical point spread function, (b): 
Regular imaging shows a ring structure, but the individ- 
ual emitters are not resolved, (c): the same image as in 
(b), smoothed by a Gaussian filter to reduce the noise (filter 
FWHM is shown by the scalebar). (d): the second order an- 
tibunching image, (e) the same image smoothed using the 
same filter as in (c). The individual emitters are clearly dis- 
cerned. The image was formed from 1000 frames. 



The antibunching signal of order n can then be defined 



as 



A r . 



(8) 



which can be expressed via cumulants of order k < n 
using Eq. ©: 



A n {x) 



k 



(9) 



where (R n )ik is the first row of the n-th power of a matrix 
R, all elements of which are zero except 



R. 
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l/k. 



(10) 



The signals A n (x) vanish in the classical limit, and are 
therefore a valid local measure of the degree of antibunch- 
ing. For n = 2 the above expressions yield the second 
order antibunching signal described above. Substituting 
eq. (fTJ) into the definition (J5J one obtains 



A n (x) = MQ n S n Y,Pl - Xi )] 



(11) 



Effectively, the n-th order antibunching signal corre- 
sponds to a point spread function h m (x) = [h(x)] m . Sim- 
ilarly to the second order, this enables imaging with res- 
olution up to n times better than diffraction limited in 
three dimensions pi |28| . 
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Quantum fluorescence imaging requires fluorophores 
with as high quantum yield as possible. Fortunately, 
many of the fluorescent markers widely used in bio- 
imaging, such as organic dyes, have a quantum yield 
approaching unity. Very high quantum yield has also 
been demonstrated with colloidal semiconductor quan- 
tum dots [Hj|. 

Many fluorescent single photon emitters exhibit ran- 
dom variations of brightness known as blinking. Blink- 
ing increases the observed photon number fluctuations, 
and could be expected to change the photon statistics 
to super-Poissonian. It turns out, however, not to be 
the case: indeed, for a single emitter, the photon num- 
ber follows a Bernoulli distribution even in the presence 
of blinking, the only consequence of which is an effec- 
tive reduction of the emission probability. As long as 
blinking of individual fluorophores is uncorrelated, the 
antibunching properties of emission statistics persist for 
an arbitrary number of emitters. 

The photon number distributions required for comput- 
ing the cumulants in Eq. ([9]) can be determined with a 
scanning number resolving detector 3^, 31 1 . Another op- 
tion is using a regular or an electron multiplying charged 
coupled device (ccd) [l5| • An electron multiplying ccd in 
the photon counting regime can be used as a pixelated 
number resolving detector, provided that the number of 
pixels sufficiently large so that the probability of detect- 
ing more than one photon in a pixel is small. Interest- 
ingly, a regular ccd, which is naturally a photon num- 
ber resolving detector, can also be used despite signifant 
amount of noise in the output. Since the noise is statisti- 
cally independent from the number of photons, it can be 
removed from the signal by subtracting the pre-measured 
cumulants of the noise from the observed cumulants. 

The antibunching microscopy can be regarded as a 
quantum version of superresolution optical fluctuation 
imaging (SOFI)[7|: in this technique, the n-th order sig- 
nal is given by n-th cumulant, without the lower order 
'correction' terms appearing in Eq. Interestingly, 
the antibunching signal vanishes in the classical limit, in- 
stead of turning into the corresponding SOFI expression. 
This is the case because the two schemes exploit different 
sources of non-Poissonian statistics: while SOFI quanti- 
fies the super-Poissonian brightness fluctuations of essen- 
tially classical sources, in the present scheme the signal 
is due to the reduction of the quantum fluctuations with 
respect to the classical shot noise. The antibunching sig- 
nal is thus generated by steadily emitting fluorophores, 
which enables continuous superresolved monitoring of the 
samples stained with fluorescent markers. 

In conclusion, we propose a fluorescence microscopy 
imaging modality that allows for sub-diffraction-limited 
imaging by virtue of quantum properties of fluorescence 



emission. Despite being ostensibly quantum, the tech- 
nique does not require a non-classical light source and 
does not depend on fragile quantum interference effects. 
The proposed method can be implemented with current 
technology, or indeed with a regular fluorescence micro- 
scope. 
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